clear;clc;close all;
figure(1)
x = 0:pi/100:2*pi;
subplot(2,2,1)
y1 = sin(x);
plot(x,y1,'r')
xlabel('x');ylabel('y1=sin(x)');
grid on
subplot(2,2,2)
y2 = cos(x);
plot(x,y1,'b')
xlabel('x');ylabel('y2=cos(x)');
grid on
subplot(2,2,3)
y3 = sin(x).*cos(x);
plot(x,y1,'g')
xlabel('x');ylabel('y3=sin(x).*cos(x)');
grid on
subplot(2,2,4)
y4 = sin(x)+cos(x);
plot(x,y1,'m')
xlabel('x');ylabel('y4=sin(x)+cos(x)');
grid on



figure(2)
x = -2:0.1:2;
y = -2:0.1:2;
[X,Y] = meshgrid(x,y);
Z = X.^2.*exp(-(X.^2+Y.^2));
mesh(Z)
xlabel('x');ylabel('y');zlabel('z');




sum = 1;
for i = 1:1:20
    sum = sum*i;
end
sum
